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The Typical Medium Theory provides conceptually the simplest order parameter description of 
Anderson localization by self-consistently calculating the geometrically-averaged (typical) local den¬ 
sity of states (LDOS). Here we show how spatial correlations can also be captured within such a 
self-consistent theory, by utilizing the standard Landau method of allowing for (slow) spatial fluctu¬ 
ations of the order parameter, and performing an appropriate gradient expansion. Our theoretical 
results provide insight into recent STM experiments, which were used to visualize the spatially 
fluctuating electronic wave functions near the metal insulator transition in Ga\- x Mn x As. We show 
that, within our theory, all features of the experiment can be accounted for by considering a model of 
disorder renormalized by long-range Coulomb interactions. This includes the pseudogap formation, 
the C(R) ~ T form of the LDOS autocorrelations function, and the C ~ i; energy dependence of 
the correlation length at criticality. 


Introduction- Until very recently, many experiments on 
metal-insulator transition concentrate on describing the 
bulk behavior averaged over the entire sample instead of 
local features on the atomic scale. That is why, the most 
theoretical approaches so far focused on global quanti¬ 
ties, such as the sample-averaged transport and thermo¬ 
dynamic behaviors, which are insufficient to properly in¬ 
terpret the new experimental data. However, a new gen¬ 
eration of STM experiments are suggesting a deep under¬ 
standing of the features of the local electronic wave func¬ 
tions in the vicinity of the metal-insulator transition^. 
Visualizing the puzzling features of the experiment on 
the nano-scale encouraged us to seek for the theoretical 
framework to describe the nature of metal insulator tran¬ 
sition. For example, one of the spectacular unexpected 
result on Gai_ x Mn x As was the emergence of pseudo¬ 
gap at Fermi energy, at the critical concentration, that 
defy the commonly accepted picture of how disorder af¬ 
fects the electronic states in the critical regiorP. This 
feature of local quantum wave functions calls for a com¬ 
pletely new look at the mechanisms for localization in the 
presence of interaction. The challenging question here 
is what type of theoretical framework can be sufficient 
and effective in order to describe the experimental fea¬ 
tures in relevant regimes. As it can be discovered by the 
experimental data, the size of the correlation length is 
comparable to microscopic cut off (the size of individ¬ 
ual manganese atoms). Moreover, the range of energies, 
which can be explored by STM measurement, is not much 
lower than the Fermi energy. Thus, the Typical Medium 
Theory such as any other mean field approaches is ap¬ 
plicable to understand the critical characteristic of this 
system^. This theoretical approach can be regarded as 
an appropriate order parameter theory for Anderson lo¬ 
calization and has the potential to directly combine with 
already well accepted Dynamical Mean-Field Theories for 
strong correlations®!!. This local order parameter theory 
has already met with considerable success in interpret¬ 
ing several experimental observation such as opening the 
pseudo-gapP, but so far were not capable to capture spa¬ 
tial correlation effects. To overcome this obstacle, we 


make a key conceptual observation guided by experience 
in standard critical phenomena. This standard method 
has been utilized for describing the variation of order pa¬ 
rameter near the critical point in the spirit of many mean 
field theories. For example, it provides a systematic way 
to understand how magnetization varies in space in re¬ 
sponse to an external field in weiss mean field theory'. 
In this letter, we construct Landau-Ginzburg theory of 
TMT for the Anderson localization by obtaining corre¬ 
sponding mean-field description of spatial correlation of 
the order parameter near the critical point. 
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Figure 1. Normalized distribution of typical density of states 
(TDOS) for uj = 0.05, 0.1, 0.2. As the transition is approached, 
the most probable value of P(p) is reduced . 

Model- Here, we study a single band tight binding 
model of non-interacting electrons with random site ener¬ 
gies ei with a given distribution P(e,), which the Hamil¬ 
tonian of this system can be written as: 

H= tijc\crCjcr + £jc\ a Cj a . ( 1 ) 

(ij),a i,<r 

Here, cj CT and Ci a are the electron creation and annihi¬ 
lation operators, and tij are the inter-site hopping ele¬ 
ments. As it has been discussed in our previous work the 
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effective disorder potential seen by quasi-particles can be 
significantly modified in presence of long-range Coulomb 
interactions . This fact leads to the formation of the soft 
“Coulomb pseudo gap” at the Fermi energy and allows 
us to consider a model distribution of random site ener¬ 
gies which assumes a pseudo-gap form expected from the 
Efros and Shklovskii (ES) pictur c P ^ l It is well known 
that the statistic of local density of states LDOS reflects 
the degree of localization of quantum wave functional, 
and its most probable, typical value, (TDOS) is repre¬ 
sented by the geometric average of LDOS. Within local 
TMT approach so far has been proven that TDOS can be 
introduced as an appropriate uniform order parameter. 
The measurement of this quantity provides the qualita¬ 
tive feature of density of states* which is in agreement 
with both experiment 1 and recent large-scale exact diag- 
onalization results 12 . As it has been shown in Fig. (1), 
near the metal insulator transition the distributions of 
(TDOS) begin to cross over from Gaussian to log normal 
distribution. Our result, which is consistent with exper¬ 
imental visualization of local quantum wave functions, 
can be another reason for justification of our simple the¬ 
ory. The missing key point in our previous work was 
ignoring the spatial correlation effects which has to be 
treated beyond the local TMT approach. The challeng¬ 
ing question is how we can constrict the nonlocal order 
parameter for Anderson localization starting with the al¬ 
ready exciting TMT we earlier proposed and tested for 
different system^ESHU] This is standardly achieved by 
introducing a weak, spatially varying external field, and 
performing an appropriate gradient expansion for the or¬ 
der parameter, which is the subject of this work and can 
be discussed as follows. 
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Figure 2. The long distance behavior of the Green’s function 
at Fermi level and also the nearby energies, r(ui) = 0.1, r(u>) = 
0.3,r(oi) = 0.5. The inset shows the same behavior in semi- 
logarithmic scale. 


propriately averaged “ effective medium” A,:(w), which 
is characterized by a self energy £(w). The self energy 
corresponds to this Landau type of TMT is local but it 
is not the same for all lattice sites. This property of self¬ 
energy is different from th e loc al TMT, which has been 
discussed extensively so f'a A 3 * 1 3 l 


cj1-H 0 -S(cj) 


I “I 


( 2 ) 


Here, H 0 is a bare lattice Hamiltonian matrix and the 
“effective medium Green’s function” Gf t M (w) is given by 
Eq. Q, which it is determined by the diagonal elements 
of lattice Green’s function. Therefore, in contrast to local 
TMT, the effective medium varies slowly in space. As¬ 
suming that the space dependent self-energy is provided 
at each site £i(w), one can calculate the “nonuniform 
cavity” by Eq.0 as 


A =w + iri- EiM ~ nEM , , • (3) 

Given the cavity field, the corresponding local density 
of states pi(uj,ei) is given by the imaginary part of the 
local Green’s function: 

Pi(u, e») = --Im[u + irj - e» - Aj(o;)] _1 . (4) 

7r 

Within Landau Ginzberg describtion of TMT, the typ¬ 
ical value of the local density of states p* yp (w) still can 
be selected as an appropriate order parameter such as 
standard TMT approach. However, the modified order 
parameter has smooth spatial variations and it is well- 
represented by the geometric average of local density of 
states as 


= exp 


daP(ei) In 


( 5 ) 


Here, in order to obey causality, the corres ponding 
“typical local” Green’s function, is defined 1 ' 18 by per¬ 
forming the Hilbert transform 


gT M = 


, p \ VP (<*>') 


du/ 

Too W + irj - u' 


( 6 ) 


The selfconsistent loop can be closed by setting the 
diagonal elements of the lattice Green function of the 
effective medium, Gf ; M (u) to be equal to the local typical 
Green’s functio n 1 H 18 !This yields to the relation 


Derivation of Landau Ginzburg theory- The general 
strategy in formulating an parameter theory follows the 
“cavity” method typically used in Dynamical Mean Field 
Theory approaches 16 . Here, the dynamics of an electron 
on a given site can be obtained by integrating out all 
the other sites, and replacing its environment by an ap- 


= ( 7 ) 

Therefore, the self-energy can be obtained by the self- 
consistency condition Eq. 0 and Eq. 0 in an iterative 
scheme as 
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Ei(w) = u + iri- Ai(w) - t . 

Gi (u) 

This method provides us an opportunity to include 
smooth spatial variation of order parameter and calcu¬ 
late it in a selfconsistent way. This Landau description 
of TMT can be obtained for any general form of density 
of states with different random distribution to explain 
the second order phase transition. This set of TMT self- 
consistent equations can be solved numerically given the 
bare lattice Hamiltonian Hq and the form of random po¬ 
tential. In this letter, we are interested in inventing the 
analytical approach to describe the critical feature of the 
system such as the behavior of correlation function and 
correlation length. The full analytical solution of TMT 
equations for uniform order parameter has been provided 
for any general modeP. Although we have found the log¬ 
arithmic correction to leading critical behavior of uniform 
order parameter at the critical point, It has been proven 
that this corrections are not important close to the metal 
insulator transition^. Knowing this fact, we are allowed 
to ignore these mild corrections to develop the non uni¬ 
form description of TMT theory. The fact that the or¬ 
der parameter varies slowly in space as the transition is 
approached, allows us to perform the gradient expansion 
for order parameter p* yp (oj). Here, we obtain the Landau 
Ginzberg equation which is given by Eq. © and provides 
us the spacial dependance of order parameter in contin¬ 
uum limit. For the sake of simplicity in notation, we set 


these has been provided for a general model of density 
of states with two classes of random distribution (generic 
and pseudogap )P Interestingly, both model provide the 
same long distance behavior for autocorrelation function. 
Within this framework, we predict that this function falls 
of! exponentially in both insulating side and metallic side 
while it shows the power law behavior at the critical point 
(W = W c ). 

For the purpose of this work, we limit out attention 
on pseudogap model of random site energies which is rel¬ 
evant to experimental data. Given this class of random 
distribution and predict all feature of experiment by solv¬ 
ing the Eq.([ 8 | exactly at the critical point and also close 
to the metal insulator transition. 
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Figure 3. The behavior of correlation length shows the diver¬ 
gence as f ~ f at critical point W = W c . 



- V 2 P + r(oj)p + Ui(uj)p 2 + u 2 (w)p 3 + ... = 0. ( 8 ) 

Here, 

r(cu) = (W-W c ) a + /3u ; 2 , (9) 

where the critical exponent a is characterized by the 
form of renormalized disorder potential. It has been com¬ 
puted for two distinct classes of random distribution as 
a = 1 (generic model, where P(w) is finite for all uj) and 
a = \ (pseudogap model, where P(w) vanishes at zero 
frequency). The value of f3 is finite at the transition. As 
it has been discussed in our previous works^DS 

ui(w) oc (P(w) + yw 2 ), (10) 

where 7 is finite and can be calculated close to the tran¬ 
sition for a given desired distribution. The coefficient 
u 2 (w) always remains finite for any model of distribu¬ 
tion. As an example, rti(u;) is finite for the generic model 
at the transition (P(0) 7 ^ 0). However, for the pseudo¬ 
gap model at the critical point (W = W c ), it is zero at 
the center of the band (w = 0 ) while it remains finite at 
nearby energies. These Landau coefficients in general are 
complicated which their form depends on both model of 
density and disorder distribution. The explicit form of 


Solving the differential equation for pseudogap model.- 
We provide the long distance behavior of autocorrelation 
function by providing boundary condition and solving 
the corresponding differential equation. We propose that 
the boundary condition is provided by having two leads 
at both metallic and insulating side. This condition is 
equivalent with exposing the order parameter to external 
field as the standard Landau procedure. The solution can 
be obtained for any model of random distribution in dif¬ 
ferent dimension. However, we focus on the model which 
is relevant to experiment and can describe the pseudo¬ 
gap formation of order parameter as the metal-insulator 
transition is approached for this model. Assuming that 
the system is uniform in three dimension, the Eq.([ 8 ]) can 
be solved for this model. Within LG description of TMT, 
we found that at finite r(w) the asymptotic behavior of 
local density of states autocorrelation function, which is 
correspond to Green’s function, shows exponential decays 
as 


G(R) 


exp{y/-rM-R} 

R 


( 11 ) 


However, at critical disorder W = W c and to = 0, the be¬ 
havior of the Green’s function is power law and is written 
as 
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G <*> = j^ = s- < 12 > 

The behavior of the autocorrelation function has been 
shown in Fig. 2. The authors also computed the angle- 
averaged autocorrelation function between two points 
separated by r. This function decays in space as 
x(Ef,r) = r~ v following power-law with r/ = 1.2 ± 0.3 
for the sample with 1.5% of Mn concentration, while, it 
falls off exponentially at all nearby energie^. Within LG 
description of TMT, we found the divergence of corre¬ 
lation length at critical point W = W c as the signature 
of metal-insulator transition. This characteristic can be 
obtained by the expression of r(w) which is function of 
distance from transition and has been shown in Fig. (3). 
This picture also it in agreement with experimental mea¬ 
surement of correlation length at critical concentration. 
The size of correlation length is in the order of two or 
three individual Mn atoms (20 A 0 ) which justifies the im¬ 
plementing of the TMT as a sufficient theory to capture 
nature of the phase transition. 

Conclusion and outlook. - In this letter, we combine the 


simplest theoretical framework, describing the disorder- 
driven metal insulator transition, with standard Landau 
procedure of order parameter. Within our invented the¬ 
ory, we predicted the several puzzling features of experi¬ 
ment on Ga\- x Mn x As. We found the pseudogap forma¬ 
tion of local density of states at the Fermi energy and 
vicinity of transition. We discovered that the LDOS au¬ 
tocorrelation functions show power-law fashion at Fermi 
energy in large distances. In addition, we verified the 
nature of second order phase transition by examine the 
corresponding correlation length at criticality. We are 
suggesting that one should measure the frequency depen¬ 
dence of entire distribution of LDOS as it can be com¬ 
puted by our theory. We claim that the entire family of 
curves which can be obtained based on existing experi¬ 
mental data, will collapse by taking our simple theory. It 
would be also thoughtful to collect experimental data on 
insulating side, where the quantum effects do not play 
very important rules, in order to be compared with our 
theoretical predictions. 
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